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1 Introduction 

PLh ' The theory of measurements continuous in time in quantum mechanics (quantum con- 

tinual measurements) has been formulated by using the notions of instrument and 
positive operator valued measure jl|]-|| arisen inside the operational approach [l], [?J to 
quantum mechanics, by using functional integrals M, M, [10|, by using quantum stochas- 



tic differential equations 11]-[23| and by using classical stochastic differential equations 
(SDE's) p^| -|Pl. Various types of SDE's are involved, and precisely linear and non linear 
' equations for vectors in Hilbert spaces and for trace-class operators. All such equations 

. contain either a diffusive part, or a jump one, or both, 

^vq In Section ^ we introduce a class of linear SDE's for trace-class operators, relevant 

to the theory of continual measurements, and we recall how such SDE's are related 
to instruments Jl], |3(| and master equations [^] and, so, to the general formulation 
of quantum mechanics. In this paper we do not present the Hilbert space formulation 
of such SDE's and we make some mathematical simplifications: no time dependence 
is introduced into the coefficients and only bounded operators on the Hilbert space of 
the quantum system are considered; for cases with time dependence see for instance 
[28, 31, ^] and for examples involving unbounded operators see for instance [p!^]-|32]. 
In Section || we introduce the notion of a posteriori state ||[ and the non linear SDE 
^ ■ satisfied by such states; then, we give conditions from which such equation is assured 

to preserve pure states and to send any mixed state into a pure one for large times. 
Finally in Section || we review the known results about the existence and uniqueness of 



invariant measures in the purely diffusive case |4C, 33, 35] and we give some concrete 
examples of physical systems. Other asymptotic results are given in 41 1. In the 
whole presentation we try to underline the open problems. 



2 Linear SDE's and instruments 

Let us denote by £(.4i;.42) the space of bounded linear operators from the Banach 
space A\ to the Banach space A2 and let us set 2(A\) = £{A\;Ai). Let H be a 
separable complex Hilbert space and let us consider the spaces of operators 2(H), 
6(H) = {pe 2(H) : Tr {p*p} < 00} (Hilbert-Schmidt operators), 1(H) = {p G £(H) : 
Tr {^/p*p} < 00} (trace-class), with the norms Halloo (operator norm) for a £ £(W), 
||p|| 2 = y / Tt{p*p} for p 6 6(H), \\p\\i = Tr{y^V} for p G 1(H). We set (a,p) = 
Tr{a*p}, for p G 1(H), a G 2(H) or for p, a G 6(H). Let us recall that 2(H) is the 
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topological dual of 1(H) and that || • ||oo < || • H2 < || • ||i- Finally, @ a quantum 
mechanical state is a trace-class operator such that p = p* , p > 0, 1 1 yo 1 1 1 ^ Tr{/o} = 1; 
we denote by S(H) the set of all states on H. We denote by || • ||i^i the norm of a 
bounded operator on 1(H) and we use similar notations for operators on the other 
spaces. 

Let H, Lj, Sh, j,h = 1,2,... , be bounded operators on H such that H = H* , 
Sfci LjLj = D\ and YlhLi ^h^h = D3 are strongly convergent in £(H). Let y 
be a locally compact Hausdorff space with a topology with a countable basis, B(y) 
be its Borel c-algebra and v be a finite measure on (y,B(y)). Let the operator 
J G £(l(H);L 1 (y,u;1(H))) be such that the maps on 1(H) defined by K A [p] 



f A J[p](y)u(dy), A G B(y), are completely positive §7|; we set = D 2 G £(W). 

According to ||, Def. 2 and Theor. 2, 7£* is a quasi-instrument. Then, we introduce 
the following bounded operators on 1(H): 



3 CO 

i=0 j=l 

1 00 

C 2 [p] = Ky[p) - -(pD 2 + D 2 p), C 3 [p] = Y, ShpS* h - -{ P D 3 + D 3 p) 



h=l 

The adjoint operators of C, C% are generators of norm-continuous quantum dynamical 
semigroups |57]]- We set also 

K[p] = C [p] + CM - l -D 2 p- 1 - P D 2 + u{y)p + CM] • (2) 

Let (p., (J-'t),^, Q) be a stochastic basis satisfying the usual hypotheses. Let Wj(t) 
be continuous versions of adapted, standard, independent Wiener processes with in- 
crements independent of the past and N(dy, dt) be an adapted Poisson point process 
on J x R_(_ of intensity u(dy)dt] N is independent of the Wiener processes and with 
increments independent of the past. We assume (O, T) to be a standard Borel space, 
(Ft) to be the natural filtration of the processes W, N and T = Vt>o J~t- 

Let us now consider the following linear SDE for T(W)-valued regular right contin- 
uous (RRC) processes: 

a t = p + / K[a s -]ds + Y] / (Lja s - +a s -L*)dW j (s) 

Jo j=1 Jo ' (3) 

+ / {j[a s -](y)-a s -)N(dy,ds), 

Jyx(o,t] 

with nonrandom initial condition p G 1(H); all the integrals are defined in the weak, 
or a(l(H),£(H)), topology of 1(H). 

The problem of the existence of a solution of eq. (|3|) can be reduced to a problem for 
SDE's in H. It is possible, in many ways, to find a larger probability space, where more 
Wiener and Poisson processes live, and to construct a linear SDE for a process ipt in H 
for which existence and uniqueness of the solution can be proved by standard means 
and such that the process at, defined by (a, at) = E[(^t|a^)] Va G £(H), satisfies eq. 
(H|) |2^, 31 1 . From this representation one obtains that there exists a solution of eq. @ 



such that: 
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i. the map p i— ► at(co) is completely positive; 



ii. if p € S(TL), then ||ot||i = Tr{a t } is a positive martingale with Eq = 1- 

About the uniqueness of the solution we are able to give some results under addi- 
tional assumptions (Theorem 1). It seems that simple estimates fail to give uniqueness 
in the general case, so this remains an open problem. Let us stress that even when 
existence and uniqueness can be proved directly from eq. (||), it seems difficult to prove 
the positivity property i. without going through the representation discussed above. 
Let us recall that the solution must be an RRC process and that uniqueness is with 
respect to Q. 

Theorem 1 (a) If Lj = 0, Vj, all the integrals in eq. ffl) are meaningful in the norm 
topology of ^(Ti), the solution is unique and also its existence follows directly from 

(V ifZr=i\\ L j\\™<+™' EZJRhWl < +oo, J[p](y) = Zn=iJn(y)pJn(yy, 

J n {y) £ £(H), f y (T^Lx \\Jn(y) WZo) v( d y) < al1 the integrals in eq. are 

meaningful in the norm topology of&(7i), the solution is unique and also its existence 
follows directly from seen as an equation in the Hilbert space &{TL). 

Proof. Case (a). All the statements follow, by standard arguments, from the estimates: 

JC[a s -]ds <\\1C\\x->i / ||cr s -||ids, 
1 ./n 



E 



/ (j[a s -](y) - a s -)N(dy,ds) < (2\\Ky\\i->i + v(y)) / E[||cr s -||i] da . 
lyx(o,t] 1J 



Case (b). The hypotheses on the operators Lj, give that K. can be seen as a 
bounded linear operator on <5(Tl); so, we have 



E 



K[a s -]ds J <t||/C|| 2 l 2 / E[\\a s -\\i]ds. 



The hypotheses on the Lj give also 



E 



[|| Ejf to*-- + ^- L j) dW ^ s )\\2 ^ 4 E n L iii- J 



E[||^|| 2 2 ]d s 



and the hypotheses on J give 



E 



yx(o,t] 



(j[a s -}(y)-a s -)N(dy,ds)\ 



< 2v{y)t 



u(dy) + uQ>) 



E Ik 



s~ 112 



Ids. 



Then, the uniqueness and the other properties follow by standard arguments. □ 

Let us denote by of the solution of eq. (||) with initial condition p (if uniqueness does 
not hold, by of we mean precisely that solution which has been constructed through a 



(4) 



SDE on H). Then g || the equation 

lt(F)[p] = Kq[1 f a?} VF € T t , Vp € S(H) , 
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defines a (completely positive) instrument [36] X% with value space (0, F t )\ the expecta- 
tion of a trace-class operator is defined in the weak topology of T(W), as the integrals in 
(|3|). An instrument X is a measure such that X(X) is a completely positive map on T('H), 
(a,I(-)[p\) is cr-additive and (l,I(Q)[p\) = (1, p); often it is the adjoint map of X which 

Instruments give the most general setting 



is called an instrument, as in [36, [28|, [HJ]. 
for representing measurement procedures in quantum mechanics: (l,I(F)[p\) is the 
probability of the event F given the premeasurement state p and X(X)[p]/||X(X)[p]||i 
is the postmeasurement state, conditional on the event F. In our case we have a time 
dependent family of instruments, but by property ii. the resulting probability measures 
are consistent and we can define a unique probability P p on (^,X") by 



P p (F) = Tr{l t (F)[p]}=E Q [K||i If] 



VXe T t 



(5) 



The interpretation of eqs. (f|) and (0) is that {Xj, t > 0} is the family of instruments 
describing the continual measurement, the processes W, N represent the output of this 
measurement and P p is the physical probability law of the output. Something more 
could be said on the instruments It, linked to some Markov property of g% (see [31] 
Proposition 2.1). Moreover, if we set 



m 3 (t) = (L J+ L*,p t -), X(t;y) = (1, J[p t -](y)) , 
it turns out |28|, [3l]] that under the physical law P p the processes 



Wj(t) = W, 



r j(t) - [ mj 
Jo 



(s)ds 



(6) 



(7) 



are independent standard Wiener processes and N(dy, dt) is a counting point process 
with stochastic intensity X(t; y)v(dy)dt. 
From eq. (JsT) it follows that 



(8) 



is the state to be attributed to the system at time t if the output of the measurement 
is not taken into account or not known; it can be called the a priori state at time t. It 
turns out that the a priori states satisfy the master equation 



d_ 
~dt 



7] t = £[r]t 



Vo = p. 



(9) 



Master equations have been introduced in quantum mechanics as a way to represent 
the dynamics of an open system, when "memory" effects are negligible. 



3 Nonlinear SDE's and a posteriori states 

If we introduce the random states 



Pt -- 

then we have, VX £ X~t, 

UF)[p]=E Q [l F o-?]=E Pp 



T P \\ ' 

h Hi 



t 111 



Ep p [iFPt] ■ 



(10) 



(11) 
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According to []39[, pt(u>) is a family of a posteriori states for the instrument It and the 
initial state p, i.e. pt{u>) is the state to be attributed to the system at time t when the 
trajectory to of the output is known, up to time t. Note that rjt = Eg [erf] = Ep ■ [pt]- 
It turns out p8| , [H|] that the a posteriori states satisfy, under the physical law P p , the 
nonlinear SDE (again the integrals are defined in the weak topology of the trace-class) 



pt = p+ I C[ps-} ds + V / (L jPs - + p s ~L* - m j (s)p s -)dW j (s) 
Jo j=1 Jo 

+ j yx (^y J[Ps-]{v) ~ Ps~) (N(dy, da) - X(s; y)v(dy)ds) , 



(12) 



where y£ = {(y,s) : y 6 3^, < s < t, A(s;y) > 0}. Let us stress that eq. ( |i~2| ) has 
a solution by construction, but that uniqueness is again an open problem, even under 
some restrictive assumption; the point is that now the law of the noises depends on the 
solution of the equation. 

An interesting problem is to look for conditions ensuring the a posteriori states pt 
for the instruments It to be almost surely (a.s.) pure when the premeasurement state 
p is pure; we recall that in the convex set S(Ti.) the pure states are the one-dimensional 
projections. An interesting information-theoretical characterization of instruments pre- 



serving pure states has been given in [42] and the structure of this class of instruments 
has been found in J43|| . A measure of "purity" of a state p is the so called linear entropy 
Tr{p(t — p)}; this quantity always belongs to the interval [0, 1) and it is if and only if 
p is a pure state. In our problem we shall consider the linear entropy g(t) = (1 — pt, pt) 
and the mean linear entropy G(t) = Kp p [g(t)] of the solution of (|i~2]). For every random 
statistical operator p t we have < G(t) < 1; moreover, G(t) = if and only if p t 
is a.s. a pure state. So the study of the behaviour of the linear entropy is a way to 



analyze whether the SDE ( |i2[ ) preserves pure states or not. As in [35 1 we can prove the 
following proposition. 



Proposition 2 Equation (li) is ensured to preserve pure states, in the sense that pt 
is a.s. a pure state for every a.s. pure initial condition, if and only if £3 = and there 
exist a set A C y, A £ B(y), a family of one- dimensional projections P y , y £ A, an 
operator V £ £(ft; L 2 (y, v; H)) such that 

(E a *«\(yu a )(y))((Vu a )(y)\, y£A, 

J[p)(y) = { (13) 

(Tr {J[p](y)}P y , y£ A, 

where p = Yl a ^a\u a ){u a \ is a spectral decomposition of p. 

Proof. We observe that, by considering the embedding of T(W) into &(Tt), the process 
pt can be viewed as an Hilbert-space valued semimartingale (cf. fl44|| , Def. 23.7). Then 
it is possible (cf. S], Theor. 27.2) to apply Ito's formula to the linear entropy g(t) and 
to the mean linear entropy G(t) as done in Proposition 1 of [0q| . In this way we get the 
conditions £3 = and J[p\{y) / Tt {J[p]{y)} is (Tr{ l 7[p](y)}z/(dy))-a.s. a pure state 



for every pure state p. Then, by using Theorem 2 of ||] and the techniques of [43], we 



get the representation (|i~3[), □ 



Let us take now equation ( |i~2] ) under the hypotheses of Proposition |2| to guarantee 
that the equation preserve pure states; in [35] we have also studied if it is possible to 
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assure also that ( |l2|) map asymptotically mixed states into pure ones. Some examples 
of this behaviour in the case of linear systems are given in []34]]. We have found some 
results only in the finite-dimensional case, the general case is still an open problem. 
The proof of the next theorem is essentially as in [pq] . 



Theorem 3 Let eq. (li) preserve pure states and let 7i be finite- dimensional. If it 
does not exist a bidimensional projection P such that 

' P{L 3 + L*)P = ZjP Vj 

< PJ*[l](y)P = q(y)P i/-a.s. (14) 
k PJ*[l](y)P = z^U-a.s. 



for some complex numbers Zj and some complex function q(y), then eq. [Vk) maps 
asymptotically, for t — > oo, mixed states into pure ones, in the sense that for every 
initial condition the linear entropy vanishes for long times: 

lim (]1 - p t , pt) = 0, a.s. (15) 

t— >oo 

The space 7i being finite-dimensional, the vanishing of the linear entropy is equivalent to 
the vanishing of the von Neumann entropy; so ( |l5|) is equivalent to lim^oo — Ti{pt In pt } 
= a.s. @. 

4 Invariant measures 

Another interesting problem is the study of existence and uniqueness of an invariant 
measure for SDE (p"2|). Here we treat the purely diffusive, pure state preserving case, 
i.e. J[p](y) = p (which implies £2 = 0) and £3 = 0; an example in the jump case 
is studied in |35|] . We assume also conditions (|14|) to hold, in order to guarantee 



that eq. (12) map asymptotically mixed states into pure ones. In this situation the 
support of any possible invariant measure is a subset of the set of the pure states 
M = \p € T(W) : p = p*, p 2 = p, Tr p = l} . Moreover, we consider only the case of a 
finite-dimensional Hilbert space; in such a case M is a compact subset of T("H). Since 
we work on a compact manifold, it is useful |45| to use SDE's of Stratonovich type. 

For dimW = n, %{TL), &(J~L) and £(H) coincide and reduce to the space of n x n 
complex matrices. By taking a suitable orthonormal basis, it is possible to represent 
every state as a (n 2 — 1) -dimensional real vector and to represent M as a 2(n — 1)— 
dimensional compact manifold imbedded into W 1 ~ 1 [35]. By working with such a 



representation of states it is possible to obtain the Stratonovich form of eq. (|12|) ; going 
back to the operator form we obtain 

00 

d Pt = A( Pt )dt + Y^B 3 { Pt ) o dWj(t), (16) 
i=i 

00 

A{ P ) = -i[H, p] + Y, {(Lj + L*, p)Bj(p) + 
3=1 

- \ [(Lj + L*)L jP - (Lj + L*,L jP )p) + (17) 



\ [ P L*(L 3 + Lj) - (Lj + Lj, P L*)p] }, 
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Bj(p) = L jP + pL* - (Lj + L*,p)p. (18) 

Equation ( |l6| ) is a Stratonovich SDE on the Coo manifold M; the A, Bj are Coo-vector 
fields. Moreover, it is possible to show that, by some rotation on the Wiener processes, 
only a finite number of Wiener processes really enter eq. (|i~6|) (this is true because 7i 
is finite-dimensional). 

In (3^] and two theorems are given about existence and uniqueness of an in- 
variant measure; both the theorems are modifications of some results in control theory 



[46]. For the notion of positive orbit of a point for the deterministic control system 



associated to a diffusion see [46] eq. (1.3). 



Theorem 4 |35| Let M be a d- dimensional Coo real compact manifold and let A, 
B\,..., B m be Coo vector fields; we consider the diffusion process defined by a Stra- 
tonovich SDE of type 

m 

dX t = A(X t )dt + Bj(X t ) o dWj(t). (19) 
i=i 

If the generator L of the diffusion process X t (L = A + \ ^jL\BjBj) is elliptic for 
every x G M\S, where S = {x\, . . . , x p } is a finite set and if in the points of S the Lie 
algebra generated by the vectors fields A,B\, . . . ,B m is full, i.e. on S 

dim C A{B 0l ... ,B m } = d, (20) 



then there exists an unique invariant probability measure p for equation frltf ) and 
Supp/u = M. 



Theorem 5 |33| Let M be a d-dimensional Coo real compact manifold and let Bq,Bi 
be Coo vector fields; we consider the diffusion process defined by a Stratonovich SDE of 
type with m = 1 . If there exists a smooth hypersurface T in M such that the field 
A is transversal to T and the equation Xt = B\{Xt) is a contraction outside V (i.e. 
3 xo € M\T such that Xt tends to xq as t — > 00 for all initial points from M\T), then, 
if we call Mq the closure of the positive orbit of xq, all the points in Mq are recurrent 
and all the points in M \ Mq are transient. Moreover there exists an unique invariant 
probability measure p and Supp/i = Mq. 

Let us observe that Theorems |] and || deal with different cases, in fact the number 
m of fields entering the diffusive part of the SDE is 1 for Theorem ||, while for Theorem 
Q one necessarily needs m > d. 

In order to check the hypotheses of Theorem || we have firstly to prove the ellipticity 
of the diffusion matrix on M but a finite set of points; the check of the ellipticity can 
be reduced to some properties of the operators Lj. Firstly, the tangent space to M in 
a point p is T p = {r S : r = r* , pr + rp = r}; moreover for p G M it is possible 

to write p = with ip G H, \\ip\\ = 1 and then it is possible to prove that we have 

also 1,^1 = {r G %(H) : r = \^){(f>\ + (f> G H, (<f>\1>) = 0}. Secondly, by using 

again the representation of the states in R n , one proves that the diffusion matrix is 
elliptic in a point p if and only if Vr G T p there exists j such that (t, Bj(p)) ^ 0. Finally 
the ellipticity hypothesis of Theorem || becomes: the diffusion matrix of the SDE (|l6|) 
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is elliptic in a point if an d only if for all £ 4> 7^ 0, (</>|V') = there exists j 

such that 



Rei^Ljip) / 0. (21) 

For what concerns the application of Theorem [| to eq. ( |l~6| ) , one has to study the 
equation p\ = B\{pt). If the initial condition is po = |V , o)(V'o|) "00 S HV'oll = 1 ; then 
the solution of such an equation can be written as 

Pt= I, ; no > % = e 1 Vo- (22) 



In [33 1 some concrete examples are studied with L\ selfadjoint and, in particular, with 
L\ one-dimensional projection. Note that eq. ( |l6|) is invariant for a change of sign of 
the Lj's and the Wiener processes; so, to check the hypotheses of Theorem [| also the 
equation p\ = —Bi(pt) can be considered. 

Even for TL = C 2 , there are some physically interesting examples. Let us consider 
a two-level atom of resonance frequency ujo, stimulated by a monochromatic laser of 
frequency u>; the interaction between the atom and the electromagnetic field is mediated 
only by the absorption/emission process. It is known that the detection schemes of the 
fluorescence light known as homodyne/heterodyne detection and direct detection can 
be treated by SDE's as eq. (|T2j), of diffusive and jump type respectively [^6[ 27, 22]. 
In the case of an homodyne/heterodyne detection scheme with a local oscillator in 
resonance with the stimulating laser, the equation for the a posteriori states turns out 
to be of the form (|l2]), with only the diffusive part and with 

Lj = (ej|a)<r_ , H = — -Acoa z + i(A|a)cr_ — i(a|A)o" + , (23) 

where the sigma's are the Pauli matrices 

1\ /() (>\ (\ 



a+ ' \o o) ' a = \i o) ' az= yo -l 

Au = to — ujq, A, a £ fC, which is a separable complex Hilbert space, {ej, j = 1, 2, . . . } 
is a c.o.n.s. in /C; the choice of this c.o.n.s. depends on the concrete implementation of 
the measurement scheme. The quantity ||a|| 2 represents the natural line- width of the 
atom and 2|(a|A)| is known as Rabi frequency; we assume both to be strictly positive, 
i.e. \\a\\ 2 >0,n = 2\(a\X}\ > 0. 

In [^] we have proved that, if the complex numbers (ej\a) are not all proportional 
to a fixed one, then eq. (|2l| ) is satisfied in the whole manifold M but in the point 

po = ( n J, where the condition (|20|) is satisfied; so, in this case, the hypotheses 



of Theorem ^| are fulfilled and there exists a unique invariant measure p which has 
Supp^ = M. 

In [ 26 1 the problem of the invariant measure for a two- level atom is studied by means 



of numerical simulations. The authors consider a heterodyne detection scheme which 
satisfy all the hypotheses of Theorem |3], as we have shown in ]35], and a homodyne 
scheme which only in a particular case can be handled via Theorem |i ; however, Theorem 
|5| can be applied. The model is characterized by m = 1, L\ = e~ 1 ^ ||a||cr_, 4> G [0, 27r), 
(a|A) = if2/2. In this case the solution of the equation p t = B\{pt) can be easily 
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computed via eq. (22) and one obtains that the hypotheses of Theorem |5| are satisfied 

with r = 0, xo = p = This means that there exists a unique invariant 

probability measure p with support Mq implicitly given in Theorem ||. From the 
structure of the fields A and B\, we conjecture that the support of p should be the whole 
M, apart from exceptional cases such as the one studied in p5| l, which corresponds to 
Au = 0, 4> = 7r /2 and for which Mq reduces to a circumference (for TL = C 2 , M can be 
represented as a spherical surface). 



Let us end by discussing the connections between the invariant measure for eq. (12) 
and the equilibrium state of the master equation (^) . The following considerations apply 
both to the diffusive and to the jump case; we assume that 7i is finite dimensional and 
that there exists a unique invariant probability measure p and that Supp p =: Mq C M. 
By definition of invariant measure we have 



W[f{p t )]p{dp) = / f(p)p(dp) (24) 
M JM 

for every bounded measurable complex function / on M; E p is the expectation in the 
case the initial condition for the process is p. Moreover, by the uniqueness of p we have 



ergodicity |48|, i.e. 

im — / 

T- 



lim i [ T f( Pt )dt= [ f(p)p(dp) a.s. (25) 

->+oo 1 J Q J Mq 



By applying these two facts to the function f(p) = (a,p), where a £ £(W), an d by 
recalling that E p [pj] = r]t satisfies the master equation (||), we obtain that the equation 
C[rj\ =0 has a unique solution solution r] eq in S(TC), given by 



Veq = / PM(dp); (26) 
JMq 



moreover for every initial condition 

,-r 

lim — / 

T 



1 f T 

hm — / p t dt = r/ eq a.s. (27) 



Another quantity of interest is D 2 (a;p) = ((a* — (a,p))(a — (a*,p)),p) = (a*a,p) — 
\(a,p}\ 2 , a G £(H), p G <S(H), which is a quantum analogue of the variance. Note that 
we have the decomposition 



D 2 (a;r leq )= D 2 (a; p)p{dp) + / \(a, p - Veq )\ 2 p(dp) (28) 

J M JMq 



and, by the ergodicity, we obtain 

f-T 

lim — / 

T 



lim 1 f T D\a- Pt )dt= f D 2 (a-p)p(dp). (29) 

^+oo 1 J Q J M 



For the case m = 1, L\ = L\, the asymptotic behaviour of D 2 {L\; p t ) has been studied 
in ||§. 
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